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The spectral properties of a disordered system with few interacting three-dimensional spinless 
fcrmions are investigated. We show the existence of a critical spacings distribution which is in- 
variant upon increase of the system size, but strongly depends on the number of particles. At the 
critical point, we report a substantial decrease of the degree of level repulsion as the number of 
particles increases indicating a decrease of nearest level correlations associated with the sparsity of 
the Hamiltonian matrix. 
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The presence and strength of a disorder potential 
strongly modify the influence of electronic interactions 
on localization. While in clean systems repulsive interac- 
tions always localize, their effect can be reversed by the 
presence of a strong disorder. There, interaction-induced 
hoppings between localized states may help electrons to 
overcome the random potential thus reducing localization 
. In this paper we will show that repulsive interac- 
tions can indeed trigger a metal-insulator transition in a 
few-particle strongly disordered system, thereby turning 
a one-body insulator into a many-body metal. 

The statistical properties of spectra of disordered one- 
particle systems are known to be closely related to the 
localization properties of the corresponding eigenstates 
[0-0]. In the localized phase, states which are close in en- 
ergy have an exponentially small overlap. Consequently 
the levels Ei are uncorrelated and the corresponding nor- 
malized spacings a* = (-Ei+i — Ei) / {Ei + \ — Ei) are dis- 
tributed according to the Poisson distribution Pp(s) — 
cxp(— s). On the other hand, the large overlap of delo- 
calized eigenstates of neighboring energy induces correla- 
tions in the spectrum and leads to level repulsion, so that 
in the time-reversal symmetric case, the distribution of 
energy spacings corresponds to the Wigner-Dyson distri- 
bution Pwd(s) = § sexp(-7rs 2 /4) M. While these two 
limiting behaviors set in progressively as the linear sys- 
tem size L increases, a third, size-independent distribu- 
tion appears at the metal-insulator transition due to the 
divergence of the correlation length j|] . This behavior 
suggested a new powerful method to locate the metal- 
insulator transition in one-particle models, relying only 
on the knowledge of the spectrum to deduce the local- 
ization properties of the cigcnfunctions. It is one of the 
purposes of the present letter to generalize this method 
to many-body interacting disordered systems. 

Spectral properties of many-body systems have already 
attracted some interest. Early works showed that, for 
clean systems, only non-generic models where special 
group symmetries ensure the integrability, have a Pois- 
sonian level spacings distribution j7j. In the presence 



of disorder, it has been suggested that the interaction 
induces correlations between many-body wave functions 
and in the corresponding spectrum, thus resulting in 
level repulsion. This phenomenon has been related to 
a derealization in real space in the Coulomb glass || 
and in a system of two interacting particles in a disor- 
dered potential [0 , and to the recently intensively inves- 
tigated delocalization in Hilbert/Fock space |0-O|. The 
latter concept refers to an interaction-induced spreading 
of the many-body states over the basis of Slater deter- 
minants built out of the one-particle states. For suffi- 
ciently large interactions, this spreading covers many ba- 
sis states whose number represents a localization length 
in Hilbert/Fock space much in the usual sense. The 
local spectral density typically acquires a Breit-Wigner 
shape pbw(E) = T/(2n(E 2 + T 2 /4)) and the spreading 
of many-body states may be characterized by the width T 
H[l2| . The latter may be estimated from the Golden Rule 
which gives T ~ Q 2 / A c where Q 2 is the mean square 
interaction matrix element and A c the spacings between 
states directly coupled by the two-body interaction. Thus 
the interaction starts to mix directly coupled states as 
the critical threshold Q c ~ A c is reached and simultane- 
ously, spectral correlations result in level repulsion and 
the emergence of quantum chaos |10 . Quite surprisingly, 
it was shown in jllj that the number of basis states oc- 
cupied by a single interacting state is given by £ ~ r/A n 
where A„ is the n-particle mean spacing. As quantum 
chaos starts to set in one has £ c ~ A c /A„ 3> 1 so that 
contrarily to the common belief, level repulsion appears 
after many noninteracting states have been mixed by the 
interaction. Since the models studied in JT^JTTJ] contained 
no information on the real space, no conclusion could be 
drawn as to a transition in real space associated to the 
emergence of quantum chaos or the mixing of eigenstates, 
a question which is of particular interest. 

In this letter, we present a systematic investigation 
of the spectral properties of a system of few interact- 
ing quasiparticles above a frozen (i.e. noninteracting) 
Fermi sea. The neglect of interaction processes involv- 



1 



ing one-particle states located in the frozen Fermi sea is 
advantageous in two respects : (i) it allows to reduce sig- 
nificantly the Hilbert space volume and consequently to 
study large systems of linear sizes up to L = 10 in three 
dimensions (Presumably, it is the impossibility of study- 
ing systems with large enough sizes that hindered inves- 
tigations similar to B in interacting systems so far.) and 
(ii) it projects out one-body tail states which are known 
to be strongly localized and hence much less sensitive to 
the interaction- induced delocalizing effect. For the two- 
quasiparticle case, this model was originally proposed in 
and further studied in |i~3| ]. 

More specifically, we consider spinless fermions in a 
disordered Anderson-like model with a repulsive nearest 
neighbor interaction of strength U whose Hamiltonian is 



H = P h 



tola* 



p b (i) 



Here, Pb projects out Slater determinants containing one- 
particle states of energy lower than the threshold energy 
Ei,. (i; j) restricts the sum to nearest-neighbor sites on 
a three-dimensional cubic lattice and Wi is the on-site 
disorder randomly distributed between — W/2 and W/2. 
In our numerical investigations we fixed Ef, «j — 3.3i cor- 
responding to a filling factor of v = 1/3 in the disorder 
regime we will consider W — 18. We point out that 
this choice is arbitrary and checked numerically that the 
physical picture presented here does not depend on this 
particular choice of E^. For a system of linear size L 
and n quasiparticles, the ground-state energy is approx- 
imately e g w nE b + n(n + l)Ai/2 + 3Un(n - l)A 1 /B 1 . 
Ai « Bi/L d is the one-particle average spacings and 
B>\ w 2dt + W is the one-particle bandwidth. We also 
define the many-body excitation energy as Se = e — e g . 
As the introduction of the frozen Fermi sea may appear 
somehow artificial, we stress that for low many-body en- 
ergy excitations Se and moderate interaction U <C B\, 
only particles in an energy layer of size T ~ t/SeAi 
around the Fermi level effectively interact and their num- 
ber is given by n e ff ~ T/Ai ~ y/ Se/ Ai so that for 
n e ff < n our model is fully justified. 

The numerical investigations were carried out much 
in the same way as in reference |5|,[l(J. We computed 
the distribution P(s) for the many-body level spacings 
Si = (e.j+i — e,)/A„ for n = 2 — 6 quasiparticles, at 
fixed excitation energies and disorder, for different in- 
teraction strengths and linear system size. As the in- 
teraction strength increases, more and more states are 
coupled and we expect quantum chaos to set in, i.e. the 
level spacings distribution undergoes a crossover from 
Poisson to Wigner-Dyson distribution. To quantitatively 
study this crossover we computed the value of the param- 
eter n = f°° [P(s) - P WD (s)]ds/ f°° [P P (s) - P WD (s)]ds. 
Here, So ~ 0.473 is the smallest root of Pp(s) = P\vd(s) 
and rj varies between 1 (P(s) = Pp{s)) and (P(s) = 



P\vd(s))- All investigations were carried out at fixed 
disorder W = 18. slightly above the critical disorder 
W c w 16.5 of the Anderson model with boxed distributed 
disorder |^,[l4|]. At this disorder, one-particle states still 
have a large enough localization length l± so that the 
short-range interaction couples a large number of states. 
We averaged over up to 10000 disorder realisations, so 
as to get at least N = 50000 levels in an energy interval 
[5e- A e /2;5e + A e /2], A £ < A x . 




FIG. 1. Evolution of the crossover parameter rj for n = 3, 
Se — 0.35, n = 4, Se = 0.5, n = 5, Se = 0.65 and n — 6, 
Se = 0.8 and disorder W = 18. (From left to right and top to 
bottom). The linear system sizes are L = 6 (circles), L = 8 
(squares) and L = 10 (diamonds). For each n, excitation 
energies satisfy n e // « 2-v/oe/Ai ~ n for L — 6. 

Fig.l shows the dependence of the crossover parame- 
ter rj as a function of U for different number of particles. 
Clearly, the situation is similar to the one encountered in 
H for the Anderson transition : the ^-curves belonging 
to different system sizes intersect at a critical U c where 
rj c is size-independent. This indicates the existence of 
a scaling due to the divergence of a correlation length 
associated to the many-body eigenstates: there is a con- 
tinuous transition from an insulating phase at U < U c 
where rj increases with the system size to a metallic one 
at U > U c where the opposite behavior is observed. We 
stress that, as already observed in different strongly dis- 
ordered systems and unlike in clean ones, the pres- 
ence of a repulsive interaction favors derealization. 

The existence of this transition dictates the value of rj c 
which consequently depends on the number of particles. 
Thus the critical distribution, i.e. the degree of level re- 
pulsion also depends on the number of particles, as shown 
in Fig. 2. r] c (n) is a monotonously increasing function of 
n and it is tempting to say that as n — > oo, the critical 
distribution tends to Pp(s), i.e. there is no level repul- 
sion at the metal-insulator transition in a system with an 
infinite number of interacting quasiparticles. 

The increase of rj c (n) with n and the corresponding 



2 



decrease in the degree of correlations between nearest 
levels may indicate that this metal-insulator transition 
proceeds much in the same way as the transition to chaos 
reported in ]l0| , i.e. that many-body derealization re- 
sults from the mixing of directly coupled states and that 
transport in the metallic phase proceeds by interaction- 
induced jumps between these states. In this case, the 
critical interaction satisfies the condition Q c ~ A c . Ac- 
cordingly, only states separated by an energy interval 
A c » A„ are correlated and the resulting level repul- 
sion can be distributed over N n = A c /A„ successive level 
spacings. Hence if N n — > oo, as n — > oo, the level repul- 
sion vanishes at the transition. 



creasing number of many-body levels, and thus explains 
how the level repulsion is reduced by an increasing num- 
ber of interacting particles. Equation (3) follows from 
the increasing sparsity of the Hamiltonian matrix with 
increasing number of particles. 

Better estimates of N n and A c should take into account 
the fact that, close to the Anderson transition, the struc- 
ture of the one-particle eigenstates reduces the number 
of effectively coupled states, while enhancing their cor- 
relations and hence the mean transition matrix elements 
p6| . The latter in their turn, are directly related to the 
spacing between directly coupled states by the relation 
obeyed at the transition Q c ~ A c . 
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FIG. 2. Critical distributions for n = 1 (squares), n — 2 
(full diamonds), n = 3 (triangles) and n = 4 (full circles). 
Also shown are the Wigner-Dyson and the Poisson distribu- 
tion (full lines). Inset : Evolution of the critical crossover 
parameter r/ c as a function of the number of particles. 

We can estimate a lower boundary for N n following a 
similar argument as in Assuming that at the tran- 
sition, the system is thermalized, we consider n — > oo 
particles at a temperature 1 <C T/Ai <C n. The number 
of effectively interacting particles is then n e // ~ T/Ai 
so that they may be distributed over m w 2n e f f different 
one-body levels. The two-body interaction connects each 
state to roughly K w u e //(n e // — l)(m— n e //)(m— n e // — 
l)/4 ~ rigj^(T/Ai) states. All these possible transitions 
occur in an energy interval B2 ~ (2m — 4)Ai « 4n e //Ai 
so that in agreement with (jnj 

B2 

K ~ nl ff T 



A, 



(2) 



Next we have A„ w 2A?(n - ly./Se 11 ' 1 (as before Se 
T 2 /Ai) and consequently N n is exponentially large 



ln(iV„) w n eff ln(n e //) 



(3) 




This latter result means in particular that at the tran- 
sition, level correlations are distributed over a fast in- 



u/u 

FIG. 3. Schematic evolution of 77 as a function o f the re nor- 
malized interaction U /U* in the regime n e ff ~ 2->J 5e/Ai < n 
(U* corresponds to the smaller system size.). The effective 
number of interacting particles increases with the system size, 
and in the process, r] c {n £ jf) increases as indicated by the ar- 
rows and the curves no longer intersect at one point. In the 
other regime studied in this article, Se > Ain, the number of 
effectively interacting particle is invariant upon increase of the 
system size and the transition is unambiguously determined 
by the intersection of all ^-curves (see Fig.l). 

Our interpretation relates the level repulsion at the 
transition to the n-body correlation volume £?i/A c , i.e. 
the volume inside which n-body states have a strong 
overlap and are therefore directly connected by the in- 
teraction. Indeed in the three-dimensional one particle 
Anderson model, this volume goes like L d_UJ with a mul- 
tifractal exponent to w 1.3 [|l7|]. Hence we have N± ~ L u 
and therefore the critical statistics is intermediate be- 
tween Poisson (JVi ~ L d ) and Wigner-Dyson [N\ w 1). 
Moreover in the case of two particles, the two-body inter- 
action connects all states together and A^2 can be reduced 
only by the one-particle correlations. Consequently, the 
critical statistics is very well fitted by the semi-Poisson 
Psp(s) = 4s exp(— 2s), as can be the case for the one- 
particle model [p.5| . As n increases, the Hamiltonian ma- 
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trix becomes sparser and sparser and the critical spacings 
distribution exhibits less and less level repulsion. 

In the thermodynamic (TD) limit, for such Fermi sys- 
tems, it is natural to define a temperature as n e ffT = 
T 2 /Ai ~ 6e, and as T must be constant as L — ► oo, this 
implies that T/B\ ss n e ffAi/Bi = const. Two different 
regimes are then of interest: 

• The one-particle critical regime close to the Anderson 
transition, where the system size is significantly smaller 
than the one-particle localization length. There, both 
the one-particle correlation volume and n e f / increase as 
~ L d ~ u until Lk.Ii and since l\ can be very large, so is 
n e ff and there is almost no level repulsion at the tran- 
sition in the TD limit. In this case the scaling at the 
transition is not observable since r\ c increases with n and 
the 77-curves intersect at different points (see Fig. 3). 

• The one-particle localized regime 1 -C h < L, where we 
have Ai as B\/lf = const. This regime corresponds to 
the results presented in Fig.l for which n e ff is constant 
and one has a fixed-point at a well-defined value of U c . 

Since it is expected that a large interaction will ulti- 
mately transform the system back into an insulator, it 
is important to estimate the parametrical width of the 
metallic phase. In this case, it seems natural to carry 
out the substitution Ai — > Q in (2) so that in the local- 
ized regime with n e tf = const, the condition Q c ~ A c 
leads to a metallic phase in the range 



Ai 



< — < 
T - Ai ~ 



<ff T 



Ai 



(4) 



Since n 2 e j^T — Se 3> Ai, the regime defined by (4) is 
parametrically large enough to allow for a metallic phase. 
The backtransition at high Q occurs before the Wigner 
crystallization, i.e. it is induced by the lowering of the 
density of states much in the same way as the transition 
described above, and occurs before the opening of a gap. 

Without the constraint of a filled Fermi sea (Pf, = 1), 
levels are filled from the bottom of the Anderson band 
where the decrease of l\ implies an increase of Ai. Cor- 
respondingly, the metallic regime determined by (4) is 
reduced unless the disorder is lowered below W c so as to 
get the one-particle mobility edge close to the Fermi level, 
and increase the localization length of the one-particle 
band-tail states. In this case, the situation should be 
similar to the one reported above, i.e. the one-particle 
correlation volume is sufficiently large to allow the mix- 
ing of a large number of states and a metal-insulator 
transition should occur. 

In conclusion we successfully applied the method devel- 
oped in H to disordered many-body systems. We showed 
the existence of a continuous metal-insulator phase tran- 
sition where the critical level spacings distribution is dif- 
ferent from the one at the Anderson transition. These 
considerations quite naturally lead to the identification 
of the metallic phase with the appearance of quantum 



chaos and the mixing of directly coupled states. 
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